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Abstract 

We consider trimaximal lepton mixing, defined by |C/a2|^ = 1/3 Va = e,/i,r. 
This corresponds to a two-parameter lepton mixing matrix U. We present a model 
for the lepton sector in which trimaximal mixing is enforced by softly broken discrete 
symmetries; one version of the model is based on the group A(27). A salient feature 
of our model is that no vacuum alignment is required. 
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1 Introduction 



It is now experimentally firmly established that the neutrinos are massive and that leptons 
of different families mix in the charged weak interaction — see [T] for reviews and [21 Bj 
for recent fits. Lepton mixing is given by a 3 x 3 unitary matrix [Uaj), with a = e,fi,T 
corresponding to the lepton flavours and j = 1,2,3 corresponding to the neutrino mass 
eigenstates. According to the most recent three-flavour neutrino oscillation update [3], at 
the 3 a level, the results for neutrino mixing are 

If/esi' < 0.056, (1) 

0.25 < sin^^Q < 0.37, (2) 
0.36 < sin^ e^tra < 0.67, (3) 

where 



1- It/, 



sin^^o = (4) 



Sin^^atm = '^l"^' 2 - (5) 

-1- — FeSi 



e3| 
|2 



Moreover, at the 3 a level the neutrino mass differences satisfy [3] 



7.05 X 10"^ eV^ < Ami = m\-m\< 8.34 x 10"^ eV^ (6) 



^0 

2 



2.07 X lO"'^ eV^ < Amft^ = - | < 2.75 x 10"^ eV^ (7) 

the sign of m| — being unknown. 

The bounds ([I])-(l3]) suggest that the lepton mixing might be tri-bimaximal, i.e. that 
the lepton mixing matrix might be, apart from unphysical rephasings and from possible 
Majorana phases, 

/ 2/V6 1/V3 \ 

t/ = t/HPs= -1/V6 1/V3 -1/v^ . (8) 
\ -1/V6 1/73 1/v^ / 

The phenomenological hypothesis U = Uups has been put forward by Harrison, Perkins 
and Scott (HPS) [Ij. However, it turns out that, at the model-building level, it is quite 
awkward to enforce U = f/nps through symmetries. In particular, all existing models for 
U = t/nps involve some form of vacuum alignment, i.e. two subsectors of the scalar sector 
having vacuum expectation values (VEVs) aligned in different directions; some papers 
where this problem has been discussed are found in [S], for more papers see the vast 
bibliography of [6J0 

In this paper we adopt the milder hypothesis that lepton mixing is trimaximal, i. e. that 

\Ue2f = K2f = \Ur2f = l- (9) 



^An alternative approach consists in using extra dimensions for model building, see for instance [7]. 
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Trimaximal mixing relaxes some of the HPS assumptions [3], since it allows for a nonzero 
Ue3 as well as for sin^ 6'atm 7^ 1/2. Our main purpose in this paper is to show that 
trimaximal lepton mixing may be enforced through a simple model which involves no 
vacuum alignment. 

In section [2] we make a brief phenomenological study of trimaximal mixing, proceeding 
in section [3] to present the simplest version of our model. In section H] we consider some 
variations on the model. Our conclusions are found in section |5l 



2 Trimaximal mixing 

It follows from the trimaximal-mixing assumption \Ue2f = 1/3 and equation (jl]) that 



1 



sin^ 9q = — > -, (10) 



which is somewhat disfavoured experimentally, since the best-fit value for sin^ 9q is 0.304 < 
1/3 \3\. The situation becomes worse for the trimaximal-mixing hypothesis when Ues is 

I 1 2 I 1 2 

nonzero; indeed, a recent fit p] found \Ue3\ = 0.016 ± 0.010, in agreement with \Ue3\ = 

0.01 t-Q^ii in [3j, which is not yet a significant indication for a nonzero Ues- In any 

case, the trimaximal-mixing hypothesis might be testable soon through more accurate 
measurements of \Ue2\ and \Ue3\. 

A trimaximal lepton mixing matrix U has the moduli of two of its matrix elements of 
the same column fixed. This means that only two parameters remain in which can 
be taken as If/gsl; or the mixing angle 613 and a Dirac phase. Clearly, for the latter a 
convention has to be adopted. Using the convention for the Dirac phase 6 promulgated 
by [H], we find 

tan2^atm = ^ ^^^^ 

\Ue3\ C0s5y 2 - 3 \Ue3\^ 

In the following, we shall employ the following parameterization of a trimaximal mixing 
matrix: 



HPS 



/ c se"*'^ \ 

1 
\ -se'^ c J 



diag{e"^\e'^\e'^') , (12) 



where c = cos 6' and s = sin 6. The mixing angle 6 parameterizes how much lepton 
mixing deviates from tri-bimaximality. The phase ip is of Dirac typejl The phases Sg^^^r, 
together with one of the phases Pj, are unphysical; only the phase differences 2 (/5i — P2) 
and 2 {P2 — P3) can be physical, if the neutrinos happen to be of Majorana type. The 
modification (fT2l) of the HPS mixing matrix has recently also been considered in |9]. 



^Besides, two Majorana phases are present in U if the neutrinos are Majorana fermions. 
■^Note that the phase tp corresponds to a Dirac phase convention different from that of 6. 
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From equations (fT^ and 



If/esi' = (13) 
. n ^ 1 CS COSIp 



.2, 1\% If/esl'l^i^ .... 
Sm ^atm - 77 < -^r- 1 \ (15) 



Therefore, 

I sin^ ^„i„ — 

V - 2 (i_|t/,3|2)' 

This inequahty can also be obtained directly from equation ( TTTl) . The inequality (IT^ 
relates, when the mixing is trimaximal, the maximal possible departure from maximal 
atmospheric-neutrino mixing, i.e. from sin^ 6'atm = 1/2, to the value of It/esl- 



3 A simple model 

Introduction Let us assume that the neutrinos are Majorana fermions. Then, in the 
weak basis in which the charged-lepton mass matrix is diagonal, the effective mass La- 
grangian for the light neutrinos is 



neutrino mass 



T T 



( VeL \ 



+ H.C. 



(16) 



where C is the Dirac-Pauli charge- conjugation matrix in Dirac space and M.^ is a 3 x 3 
symmetric matrix in flavour space. The lepton mixing matrix V diagonalizes M.^'- 



M.^ U = diag (mi, m2, m^) 



(17) 



the neutrino masses mi 2,3 being non-negative real. Using the parameterization of a tri- 
maximal U in equation (IT^ . we shall denote 



Hj = rrije for j = 1, 2, 3. 

Then, if we assume the phases 5e,^,r to vanish, we have 

/ 

diag(/ii,/i2,Ai3) = 



v 



c 
1 

se~'^ 





c 



\ 



c 


'Se 



se-''^' \ 

1 







c 



(19) 



Thus, up to the phase transformation given by the phases 5e,^t,T; trimaximal mixing means 
that the vector (1, 1, 1) is an eigenvector of A^,^ with eigenvalue fj,2- This means that, in 
the phase convention 6e = = 6r = for Aii,, the sum of the matrix elements of Ai^ 
over all rows and columns of M.^, is equal (to /i2). It is the purpose of this section to 
construct a model based on this idea. 
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S 


SU{2) X f/(l) 


(2,-1) (1,-2) 


(1,0) 


(1,0) 


(2,1) 


(2,1) 


(1,0) 


A(27) 


3 3* 


3 


1(1,0) 


3* 


1(0,0) 


1(1,0) 



Table 1: Fermion and scalar multiplets of our model 



The group A (27) A (27) is a discrete group with 27 elements. It has two inequivalent 
triplet irreducible representations (irreps), 3 and 3*, and nine inequivalent singlet irreps, 
]^(p,9) (j)^q = 0,1,2). The triplet irreps of A(27) are faithful, the singlet irreps are non- 
faithful. The group A (27) is generated by two transformations, C and T. In the 3, those 
transformations are represented as 



3 : C 



( ° 





1 ^ 




/ 


1 





\ 


1 




















Vo 


1 


oy 




v 











with uj 



,27ri/3 



■1 +iv/3 



(20) 



Notice that the matrices representing C and T belong to S'[/(3), therefore A (27) may be 
viewed as a subgroup of S'f/(3). In the 3*, 



(21) 



:22i 





/ 





1 \ 




/ 1 





\ 


3* : C ^ 


1 




















VO 


1 


J 




vo 





u J 



In the singlet irreps, 

l(p,9) : C ^ujP, T 
The relevant tensor products of irreps of A(27) are 

3* ©3* S 



3 ® 3 
3 3* 



3*, 



^2 1 (P,g) 



(23) 
(24) 



Multiplets and symmetries In our model we consider only the lepton sector and the 
electroweak interactions. The gauge group is the standard SU{2) x U{1). There are three 

rp 

left-handed-lepton doublets Dai = {^aL, c^l) and three right-handed charged- lepton 
singlets an (a = e, ^, r). We add to these standard multiplets four right-handed neutrino 
singlets in order to enable the seesaw mechanism [10]. Those four right-handed neutrinos 
are divided in two sets, three UaR and one uor. In the scalar sector, there are four Higgs 
doublets, once again divided in two sets: three ipa and one 0o. We need moreover a scalar 
gauge singlet S. The SU{2) x U{1) and A(27) assignments of all these multiplets are 
given in table [H 

Additional symmetries Besides the gauge group and A (27), we impose three extra 
Z2 symmetries '^e,^ly■ 

Zq : -ur, (pa -0a, (25) 
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while all other multiplets remain unchanged. Each Zq, has the purpose of "gluing" to 
(pa in the Yukawa couplings; this is the same idea as in [11] (see also [12j). Since the 
do not commute with A(27), the horizontal symmetry group employed in our model is 
actually much larger than A (27). 

Yukawa couplings Let us first consider the Yukawa couplings of the They are 

•^dfl Yukawas Ul 

5aL«iJ0a + H.C. (26) 

According to equation (l23l) . A(27) would allow two other couplings, 

-y2 {DeLfJ'RA + D,j.LrR4>e + DTLeR(f)^) ^^y) 
-ys {DeLTR(t>^, + D^LCRCpr + DrLfJ'RA) + H.C. 

These terms would destroy trimaximal mixing. It is for this reason that we have in- 
troduced into our model the symmetries Zq,, which remove the terms (1271) from the La- 
grangian. The masses of the charged leptons are rria = \yiVa\, where is the VEV of the 
neutral component of 0^. If we manage fe,/i,T to be all different, then the charged leptons 
will be non-degenerate in mass as desired. 

The Yukawa couplings of the right-handed neutrinos are given by 

Yukawas = "1/4 ^ D^LJ^aR («^20o) + y '^Or'^'^'^OR S + H.C. (28) 

Soft breaking of the symmetries Soft breaking of (super)symmetries plays an impor- 
tant role in many models. Soft breaking is usually needed in models which want to explain 
mixing features through some symmetries. It has been emphasized that successful models 
need a residual symmetry [13]; for instance, in [TTl IH] the residual symmetry after soft 
breaking is the fi-T interchange symmetry, which leads to maximal atmospheric-neutrino 
mixing and to Ues = 0. 

In the present model, we break A (27) softly in two steps. Firstly we allow it to 
be broken, by terms of dimension three, down to the symmetry generated by the 
transformation C, which is denoted Tj^^C). Secondly we allow Tj^^C) to be softly broken 
by terms of dimension two. We thus have the soft-breaking chain 

A(27) Z^iC) {e}, (29) 

where {e} symbolizes the trivial group consisting only of the unit element. 

The soft-breaking terms of dimension three occur in the Lagrangian of bare Majorana 
masses 

_ ^ sr^ T ^-1 ^ 

*-'Majorana masses 2 / ^ '''^aR^ '^aR 

M* 

+ ^ {/eR + ^ ^Jii + y^rR) C'^^OR + H.C. (30) 
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The soft-breaking terms of dimension two occur in the scalar potential 

V = Ce(pl4>e + C^<Pl4>fi + Cr0r0r H , 



(31) 



the coefficients (with dimension mass squared) Ce, and Cr being all different, thereby 
breaking Tj^iC). This is needed in order to obtain, upon spontaneous symmetry breaking, 
different VEVs v^. ur said, therefore, different charged-lepton masses: 



rrie : : 



Furthermore, we might include in V all terms like 



rl- (32) 

etc. This would also break the 
symmetries Zq, softly and would avoid all potential problems with spontaneous breaking 
of discrete symmetries in our model. 

Seesaw mechanism From equations (128|) and (130|) . we see that there are in our model 
Majorana and Dirac neutrino mass matrices 



R 



/ Mo 


Ml 


Ml 


M2 




Ml 


Mo 


Ml 


cj^Ma 




Ml 


Ml 


Mo 


UJM2 




\ M2 




UJM2 


Mn 


1 



M, 



D 



/ a \ 

a 

0a 

V / 



(33) 



respectively. We have defined = y^Vg, where vs is the VEV of the scalar singlet 
S, and a = ylvo, where Vq is the VEV of the neutral component of (pQ. The seesaw 
mechanism [lO] tells us that 

= -MIM^^Md. (34) 
After some algebra one finds that M.^, is of the form 



/ X + y z + uj'^y z + toy \ 



z + uj y X + ujy 



z + y 

,2„ 



(35) 



where 



X 



\ z + Ljy z + y X + UJ y J 
Mo + Mi 



—a 



(Mo - Ml) (Mo + 2Mi) 
2 ^1 



—a 



(Mo - Ml) (Mo + 2Mi) 
2 



M| 



Mn (Mo - Ml) 



2 • 



It is clear that 





( ^ ^ 




( ^ ^ 




1 


= {x + 2z) 


1 











(36) 
(37) 
(38) 

(39) 



Therefore, our model predicts trimaximal mixing. 
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A further prediction We compute 



[ X - z + \ y i|y \ 

x + 2z 

3„. n ^ ^ 3 



(40) 



Comparing this result with equation f|T9l) . we see that 



-a" 



and 



2 

Therefore, 



U2 = a; + 2;^ = — — (41) 

^ Mo + 2Mi ^ ' 



x-z + ^y = /iic^ + zissV^, (42) 
x-z-^y = /i3c2 + /iisV2*^ (43) 
z^y = cs [fise'^ - f,,e-'^) . (44) 



/ii ( c — ise^'^ 



hence 

mi 1 + 2cs sin ^p 
rris 1 — 2cs sin ip 
Comparing this result with equations ([T3l) and ([T^ . one finds that 



(45) 
(46) 



(2 - 3 - (1 - If/eaH (2 sin^^atm - l) = tt ^, (47) 

o (mi + ms) 

c/. inequahty f[T^ . The prediction fB7|) relates the deviation from tri-bimaximal lepton 
mixing to the mass ratio mi/m^. Using the experimental value of Am^^j^, then, the sum 
of neutrino masses mi + ms is determined by the deviation from tri-bimaximal mixing. 
With the experimental 3(T bounds ([1]) and ([7]), one finds 

mi +m3 > O.OeOeV. (48) 

If If/esl^ is smaller than the bound ([T]) and/or if sin^ ^atm deviates from 1/2, then the 
lower bound (HSj) on mi + ms is strengthened. However, the fourth power of mi + in 
equation fj471) dampens this effect. 

Taking the experimental values of Am^tm and Amg as input, equation P7|) determines 
the smallest neutrino mass mmin, which is mi for the normal and ms for the inverted 
spectrum, as a function of |f/e3| and sin^6'atm- In Fig. [T] we have depicted 
function of |?7e3|^, fixing sin^ 6'atm at 0.5 and using the mean values Am^^^^ = 2.4 x 10~^ 
eV^ and Amg = 7.65 x 10~^ eV^ from [3]. We also show the sum of the neutrino masses 
for both the normal and the inverted spectra. We see that at large \Uez\ the sum of the 
neutrino masses is safely below present cosmological bounds 
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Figure 1: The minimal neutrino mass and the sum of the neutrino masses, for both types 
of spectra, as a function of jf/esl • We have fixed the mass-squared differences at the mean 
values given in [3] and assumed that atmospheric mixing is maximal. 

Parameter counting and the number of predictions The neutrino mass ma- 
trix fl35l) is a five-parameter mass matrix because it has three complex parameters, with 
only their relative phases having a physical meaning. One can easily show with equa- 
tion fHOl) that AilAiu has only four parameters. The neutrino masses, the mixing angles 
and the Dirac phase follow all from AilAi^,, only for the investigation of the Majorana 
phases we need in fact J^^. 

Therefore, the four parameters in AilAi^y determine seven observables. As a conse- 
quence, there must be three predictions. Indeed, those predictions are given by equa- 
tions (fTOl) and (ITT]) , which follow from trimaximal mixing alone, and equation (H7|) . which 
is a result of our specific model. 

As for A4iy, the difference of the two Majorana phases can be expressed in terms of its 
five parameters; this constitutes the additional prediction if we consider Aii, instead of 
AilAiu- However, in our model, we expect no significant result for the effective neutrino 
mass in neutrinoless double /3-decay as compared to the general case, since one Majorana 
phase is competely free and \Ue3\ is small. 
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4 Variations on the model 



4.1 Use of a CP symmetry 



The model presented in the previous section does not possess /i-r interchange symmetry, 
which would require, in the A4i, of equation ( !35l) . y = 0, leading to two degenerate 
neutrinos. Alternatively, though, we may impose on the model a. CP symmetry which 
interchanges the /i and r families [I6], viz. 



( VeR \ 
\ ""rR ) 



( \ 



CP 



eR 

-T 
tR 



V. 



{t, -f) , uoR (t, f ) ^ 2 70 C z>o« (t, -f ) , (49) 



and so on. Such a CP symmetry would force Mo,i,2 in equation ( l30ll to be real, hence x, 
y and z in equations (!35|) - (!38|) to be real. One would thus obtain a neutrino mass matrix 
with three parameters only, which fulfils 



SMuS = Ml with S 



/ 1 \ 

1 
\0 1 o) 



(50) 



As shown in [T^, a neutrino mass matrix obeying equation fl50|) predicts sin 9, 



= 1/2, 

maximal CP violation, i.e. e^^ = ±i, and vanishing Majorana phases. Therefore, a 
restricted version of our model, including the CP symmetry fH^ . has these predictions 
in addition to those of trimaximal mixing. 



4.2 One more right-handed neutrino 

If one wants to have trimaximal mixing without the extra prediction (H7!) . one may in- 
troduce into the model one more right-handed neutrino, i/q^, transforming as under 
A (27). This leads to one extra term 



y^u',ic-x^s*+}i.c. 



on the right-hand side of equation fl28l) . and to two extra couplings 



Ms ( 



on the right-hand side of equation (1301) . Equations (!33|) would then read 



(51) 



(52) 





/ 


Mo 


Ml 


Ml 


M2 


Ms 






/ 


a 








\ 






Ml 


Mo 


Ml 


cj^Ms 


uMs 











a 







Mr = 




Ml 


Ml 


Mo 


UJM2 






, Mb = 










a 








M2 


UJ^M2 


UJM2 


Mn 


M4 





















\ 


M3 






M4 




) 




V 











/ 



(53) 
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OiR 




4>a 


T 


T 




T 


rp2 



Table 2: Generalizing T. 



with m'j^ = HqVs- Instead of equation (!35|l one would then have 

/ X + y + t z + uj'^y + ut z + uy + uH \ 



Mr 



z + uj^y + ut X + ujy + ujH 



z + y + t 

,2„ 



(54) 



y z + coy + uj t z + y + t x + to y + cut J 
This still predicts trimaximal mixing but the extra prediction (H71) disappears. 

4.3 Use of A(3n2) or other symmetry groups 

We may generalize the symmetry T by using, instead of equation (j20|) . 





(0 





1 \ 




/ 


1 








\ 


3 : C = 


1 















cr 









^0 


1 


J 




I 








a* 


/ 



with cr 



{n > 3). (55) 



The transformation properties of the multiplets under T would then be as shown in table El 
fields not shown in that table transform trivially under T. All fields transform under C in 
the same way as in section [3l Apart from T, all other details of the model are the same, 
in particular the soft breaking as given by equation fl29l) . Thus, the matrix Mr would 
remain unchanged, since its form is fixed by the transformation C. 
Let us consider n > 4. Then the main conclusions are the following: 

• The fermionic sector and, therefore, the matrix A^,^ is the same for all n, as given 
in equation (15^ . 

• The terms in equation ( 1271) are automatically forbidden, therefore the symmetries 
of equation ( l25l) are not needed. 

• The symmetry group is A(3n^), softly broken by terms of dimension three to ^3(6*). 

A detailed discussion of A(3n^) is given in [l7j. Actually, the terms of dimension four 
in the Lagrangian are invariant under all the permutations. This leads to the symmetry 
group A(6n^) — see [T7] . 

There is still another way to produce the present model. Consider cyclic permutations 
(or all permutations), plus family lepton-number symmetries U{l)a and the symmetries 
of equation (l25l) . The scalar doublets carry no lepton number. The neutrino uqr and the 
scalar singlet S may either carry lepton number or not. The U{l)a are softly broken by 
terms of dimension three, the residual symmetry being once again Z3(C). 

All the groups considered here produce identical models as far as the terms in the 
Lagrangian involving the fermion fields are concerned; the only differences which may 
arise occur in the terms of dimension four in the scalar potential. 
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5 Conclusions 



In this paper we have focused our attention on trimaximal lepton mixing, with a two- 
parameter lepton mixing matrix, which generahzes tri-bimaximal mixing. Our main mo- 
tivation was to allow for a deviation of IL^esI from zero; recent studies point out that 
possibility [21 [3]. Trimaximal mixing correlates the deviation of \Ue3f from zero with a 
deviation of sin^6'0 from 1/3 and of sin^ 6'atm from 1/2 — see equations ffTOl) and (|TT|) . A 
particular consequence is sin^ > 1/3, which is slightly disfavoured by the data at the 
moment, but in any case might be tested soon. 

We have also constructed a seesaw model (or rather a class of models) where trimax- 
imal mixing is enforced by a family symmetry group. In this model, the mass matrix 
of the light neutrinos, given by equation (!35|) . has five physical parameters; it includes 
not only the predictions of trimaximal mixing but also equation (147]) which relates the 
deviation from tri-bimaximal lepton mixing to the ratio nii/m^ of neutrino masses. As 
for a family symmetry group, we have considered several possibilities; one of the most 
straightforward ones is based on the group A (27). We have also considered a restricted 
version of our model by imposing, in addition, a non-standard CP transformation; in 
this way we are lead to a three-parameter neutrino mass matrix which predicts not only 
trimaximal mixing but also sin^6'atm = 1/2. 

Our model has some peculiarities, like the need of four (or more) right-handed neutrino 
singlets; the fourth neutrino singlet couples to the the three usual ones denoted by VaR 
(a = e, /i, r) via a scalar gauge singlet. In the mass terms of the PaRi the symmetry 
group is broken softly down to a Z3. An outstanding feature of the model is that no 
vacuum alignment is required, despite its scalar content of four Higgs doublets and the 
scalar singlet. This is to be contrasted with models for tri-bimaximal mixing which are 
plagued by the intricacies of vacuum alignment. 
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